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1  Introduction 

At  standard  temperature  and  pressure,  many  ceramic  crystals 
and  crystalline  minerals,  e.g.,  corundum  (a-AUC^),  quartz 
(o,-Si02),  moissanite  (6/f-SiC),  and  diamond,  are  hard,  brittle, 
and  elastically  stiff  and  often  exhibit  a  predilection  toward  frac¬ 
ture  over  dislocation  glide  or  deformation  twinning.  However,  at 
high  temperatures  or  high  pressures,  dislocation  glide  or  twinning 
can  be  important  [1,2],  Confining  pressures,  for  example,  those 
occurring  in  impact  loading  or  indentation,  can  partially  or  fully 
suppress  fracture  [3-5]. 

The  theory  developed  in  the  present  paper  addresses  behavior 
of  such  materials  in  regimes  wherein  nonlinear  thermoelastic 
properties  as  well  as  inelastic  deformation  modes  may  be  of  im¬ 
portance.  The  constitutive  model  developed  in  this  paper  is  di¬ 
rected  toward  cubic  Bravais  lattices  and  centrosymmetric  struc¬ 
tures.  For  example,  spinel  ceramics,  such  as  MgAl204  and 
aluminum  oxynitride  (AION),  belong  to  point  groups  with  cubic 
symmetry  and  centrosymmetry  [6],  The  general  framework  devel¬ 
oped  here  applies  to  any  crystal  with  such  symmetries  exhibiting 
the  following  physics:  hyperelasticity,  thermal  expansion,  slip, 
twinning,  and  shear  stress-induced  fracture.  The  particular  appli¬ 
cation  is  a  study  of  AION,  a  transparent  ceramic  known  to  exhibit 
dislocation  glide  and  twinning  mechanisms  when  subjected  to 
confined  compressive  loading  at  high  strain  rates  [4,7].  Disloca¬ 
tion  glide  has  been  observed  in  other  ceramics  sharing  the  spinel 
structure  (e.g.,  stoichiometric  variations  of  Mg0-«A1203,  1  Sn 
£3.5)  at  various  confining  pressures  and  temperatures  [8,9]. 

Because  single  crystals  of  AION  cannot  presently  be  obtained 
using  standard  processing  routes,  they  have  not  been  mechanically 
tested.  Anisotropic  second-order  elastic  constants  for  AION  have, 
however,  been  recently  computed  theoretically  via  first  principles. 
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i.e.,  quantum  mechanics  [10],  Polycrystalline  AION  has  been  ex¬ 
tensively  tested  and  characterized  [7,11,12],  Polycrystals  exhibit 
typical  grain  sizes  of  20-200  /rrn  with  minor  porosity  (1-5%, 
often  near-spherical  pores  of  diameter  1-15  pm)  and  impurities 
often  located  preferentially  along  grain  boundaries  [7,13].  Results 
presented  later  in  this  paper  demonstrate  how  a  model  incorporat¬ 
ing  nonlinear  thermoelasticity  and  pore  collapse  can  explain  the 
hydrostatic  (pressure-volume)  response  of  shock-loaded  polycrys¬ 
talline  AION  to  pressures  up  to  ~30  GPa  and  suggest  how  slip, 
twinning,  and  shear  fracture  may  contribute  to  the  deviatoric  (i.e., 
shear  stress-distortion)  mechanical  response.  The  model  devel¬ 
oped  here  may  be  used  later  in  mesoscale  calculations  of  defor¬ 
mation  and  failure  of  polycrystalline  microstructures  [14,15], 

This  paper  is  organized  as  follows.  Section  2  describes  a  gen¬ 
eral  constitutive  model  framework  for  cubic  ceramic  crystals  and 
minerals  exhibiting  a  center  of  symmetry.  Section  3  specializes 
the  model  to  single  crystals  of  AION  and  includes  a  compilation 
of  thermomechanical  properties  and  properties  of  lattice  defects 
such  as  full  and  partial  dislocations  and  deformation  twins.  Prop¬ 
erties  are  obtained  via  theoretical  considerations  (e.g.,  the  Peierls- 
Nabarro  model  or  anisotropic  elasticity  [16,17])  when  experimen¬ 
tal  data  are  absent.  These  computed  defect  properties  have  not 
been  reported  elsewhere.  Section  4  develops  an  averaging  scheme 
to  predict  the  response  of  polycrystals  with  random  initial  lattice 
orientations  and  initial  porosity,  incorporating  treatment  of  reduc¬ 
tion  in  elastic  coefficients  of  the  macroscopically  isotropic  aggre¬ 
gate  using  a  self-consistent  approach,  pore  collapse  at  high  pres¬ 
sures  [18,19],  and  possible  fracture  induced  by  shear  stresses  at 
grain  boundaries  or  other  defects.  Discussed  are  predictions  of  the 
response  of  polycrystals  to  hydrostatic  compression  and  of  poly¬ 
crystals  and  single  crystals  of  several  orientations  to  uniaxial 
strain  loading  under  adiabatic  conditions  at  a  high  strain  rate.  The 
latter  set  of  loading  conditions  is  representative  of  plate  impact 
experiments  on  single  crystals  and  polycrystals  [2,12,20-23]. 
Comparisons  with  experimental  data  on  polycrystals  [12,24-27] 
are  included;  results  for  single  crystals  are  limited  to  predictions 
because  single  crystal  experimental  data  are  not  available. 
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Terminology  and  notation  of  nonlinear  continuum  mechanics 
are  used.  Vectors  and  tensors  are  written  in  indicial  notation  in 
Cartesian  coordinates.  Repeated  subscripted  indices  are  summed. 
Superposed  - 1  and  •  denote  inversion  and  material  time  differen¬ 
tiation,  respectively. 


2  Nonlinear  Model  for  Cubic  Single  Crystals 

Let  xa=xa(XA,t)  denote  the  spatial  coordinates  that  depend  on 
reference  coordinates  XA  of  a  material  particle  and  time  t.  The 
deformation  gradient  FaA  is  decomposed  into  elastic  (FEa)  and 
plastic  (f/j)  parts  as  [28,29] 

FaA  =  8Axa  =  FEaaFpaA  (1) 

with  dA(. )  =  <?(.)/ dXA  =  FaA <?„(.).  Denote  by  va=xa  the  particle  ve¬ 
locity.  The  spatial  velocity  gradient  is 

Fab  =  db*a  =  dbva  =  F 'oaF-A  =  FEaaFE  lab  +  FEaLPpFE  1  pb  (2) 

The  total  plastic  velocity  gradient  results  from  dislocation  glide 
and  deformation  twinning: 

FPap  =  FPaBFP~lBl g  =  2  is'jn'n  +  2  HJjnfp  (3) 

i  i 

In  Eq.  (3),  the  set  {y,sla,m1^  denotes  {shearing  rate,  slip  direc¬ 
tion,  and  slip  plane  normal}  for  slip  system  i  and  the  set 
\JJ  ,y ,  s]a ,  mJp\  denotes  {volume  fraction  rate,  fixed  twinning  shear, 
shear  direction,  and  habit  plane  normal}  for  twin  system  j.  The 
total  volume  fraction  of  twinned  crystal  is  fT=2jf.  A  comprehen¬ 
sive  treatment  of  kinematics  of  crystal  plasticity  models  incorpo¬ 
rating  twinning  is  given  elsewhere  [2].  Since  j  and  m  are  orthogo¬ 
nal  for  each  i  or  j,  volume  is  unchanged  by  plastic  deformation: 
jp=d(det  Fp)/dt=JpLPa=0.  Local  mass  conservation  results  in 
p0=pJE,  where  p0  is  the  reference  mass  density,  p  is  the  spatial 
mass  density,  and  7£=det  FE. 

Denote  by  8  the  absolute  temperature  and  by  E^p=(FEaFEp 
-Sap)/  2  the  finite  elastic  strain  tensor  with  Sap  Kronecker’s  delta. 
Helmholtz  free  energy  per  unit  mass  ifA.EEap,8,pD,  tjt)  is  pre¬ 
scribed  as 


Table  1  Structural  features  of  y-AION 


Feature 

Value 

Remarks 

Reference 

Composition 

AI23O27N5 

35.7  mol  %  AIN 

[39] 

Crystal  structure 

Spinel 

Cubic 

[6] 

Space  group 

Fd3m 

No.  227 

Point  group 

m3m 

Centrosymmetric 

Lattice  parameter 

7.952  A 

Room  temperature 

[13] 

A=l/[e0(s-l)]  (7) 

with  e0  as  the  vacuum  permittivity  and  e  as  the  static  dielectric 
constant.  In  centrosymmetric  crystals,  piezoelectric  and  pyroelec¬ 
tric  effects  do  not  appear.  Electrostrictive  constants  unknown  for 
the  material  of  subsequent  interest  (i.e.,  AlON)  are  omitted.  The 
rightmost  two  terms  in  Eq.  (4)  account  for  energies  of  lattice 
defects.  The  scalar  dislocation  density  (length  per  unit  volume)  is 
pD  and  the  scalar  twin  boundary  density  (area  per  unit  volume)  is 
r]T.  Scalar  b  denotes  the  magnitude  of  a  typical  Burgers  vector,  K 
is  an  energy  factor  that  depends  on  elastic  constants  and  disloca¬ 
tion  type  [16],  and  y  is  the  stacking  fault  energy  presumed  twice 
the  energy  per  unit  area  associated  with  twin  boundaries  [2,17], 
The  elastic  second  Piola-Kirchhoff  stress  Sap  and  Cauchy 
stress  aab  are  obtained  from  Eq.  (4)  as 


-■JEFEcA^abF%1- 


Po 


dl/j 

8Ee 


'  F apxfiE xS  A  F (yfixbr.(i>F- xb>F-  tab 


-P(8-80)Sap  (8) 

Maxwell’s  stress  of  nonlinear  electromechanics  [33]  is  omitted 
in  Eq.  (8)  since  applications  considered  later  do  not  involve  elec¬ 
trodynamics  or  strong  electric  fields.  After  applying  thermody¬ 
namic  arguments  [2,23,31],  the  local  balance  of  energy  can  be 
expressed  as 


p0cv8  =  kX2  0  +  p0i-  +  2  i* V  +  2  Ffy1-  OfiE^ 

i  j 

sk\  .  i_ 

K  -  6 —  b~pn  -  ~  yriT 

381  2 


,E 

aa 


(9) 


Po<A=  \capx£EapEExS+  '-Fapx^pE^A- /3EEa(e-  80) 

81  -  ,  1_ 

-  p0cv8  log—  +  ~APaPa  +  Kb  pD  +  - yyT  (4) 

C/0  Z  Z 

The  second-  and  third-order  elastic  coefficients  are  Capxg  and 
Capxs£$-,  respectively;  these  coefficients  may  depend  on  tempera¬ 
ture  in  Eq.  (4).  In  cubic  crystals  of  Laue  group  Cl  [30-32],  inde¬ 
pendent  second-order  coefficients  in  Voigt’s  notation  are  Cy,  Cl2, 
and  C44  and  independent  third-order  coefficients  are  Cm,  C 112, 
C123,  C144,  C155,  and  C456.  The  isochoric  specific  heat  per  unit 
mass  is  cv  and  80  is  a  constant.  The  thermal  stress  coefficient  /3 
satisfies 


/3=3aTB  =  p„Tcv  (5) 

with  B  as  the  isothermal  bulk  modulus  (B  =  (Ci1  +  2C12)/3  at  null 
elastic  strain),  aT  as  the  coefficient  of  thermal  expansion,  and  T  as 
the  Gruneisen’s  parameter.  Denoting  by  cP  the  isobaric  specific 
heat  and  By  the  isentropic  bulk  modulus,  the  following  identity 
applies  [30]: 

Cplcv  =  BrJB=  1  +  3aTV6  (6) 

Isothermal  shear  moduli  (Cu-Ci2)/2  and  C44  are  equal  to  their 
isentropic  counterparts  in  cubic  crystals.  The  reference  electric 
polarization  vector  is  PA  =  FaApa  [33]  with  pa  as  the  spatial  polar¬ 
ization  per  unit  volume.  The  material  constant  A  satisfies  [34,35] 


with  k  as  the  thermal  conductivity,  r  as  an  external  heat  supply, 
and  the  resolved  shear  stresses 


F  =  JEFEJaaubmipF%\ 


/■OX 


J 


•bmlpFpl 


(10) 


Kinetic  equations  for  slip  and  twinning  are  prescribed  as  [2] 


y=roi///cr  sgn(A  f 

_|(ro/y)|///cl"  IV  >0  and  fT<  1] 

1  0  [/SO  or  fT=  1] 

where  y0  is  a  constant,  m  and  n  are  the  rate  sensitivities,  and  /c 
and  tJ-  are  the  slip  and  twin  resistances  that  can  evolve  with  the 
thermodynamic  state  as  well  as  activity  of  all  slip  and  twin  sys¬ 
tems.  Conditions  in  the  second  of  Eq.  (11)  require  that  the  rate  of 
twin  volume  fraction  is  non-negative  (stress-induced  detwinning 
is  not  considered)  and  that  the  total  volume  fraction  of  twinned 
crystal  /j-Sl. 


3  Aluminum  Oxynitride:  Thermoelastic  Properties 
and  Lattice  Defects 

Structural  parameters  for  AlON  with  supporting  references  are 
listed  in  Table  1 .  Thermophysical  properties  are  listed  in  Table  2; 
many  follow  from  a  review  paper  [11],  which,  in  turn,  lists  a 
number  of  sources.  Elastic  properties  are  listed  in  Table  3.  Be¬ 
cause  only  isotropic  polycrystalline  values  are  available  for 
pressure-  and  temperature-derivatives  of  elastic  coefficients  [13], 
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Table  2  Thermophysical  properties  (at  room  temperature  0O  unless  noted  otherwise) 


Parameter 

Value 

Definition 

Reference 

Po 

3714  kg/m3 

Mass  density 

[13] 

Cp 

708  +  2.32(0- 0O)  J/kgK 

Isobaric  specific  heat 

[11] 

aT 

5.23  X  10-6/K 

Thermal  expansion 

[11] 

k 

10.89  W/mK 

Thermal  conductivity 

[11] 

s 

8.5 

Dielectric  constant  (at  102  Hz) 

[H] 

Sm 

2165  K 

Melting  temperature 

[11] 

r0 

1.28 

Gruneisen  parameter 

Eq.  (5) 

such  derivatives  of  second-order  elastic  coefficients  are  estimated 
in  model  calculations  from  isotropic  results  via  CIU-C[2-2CI44, 
where  primes  denote  a  partial  derivative  with  respect  to  tempera¬ 
ture  or  pressure.  Third-order  coefficients  have  not  been  measured 
experimentally  but  can  be  estimated,  in  part,  from  pressure  deriva¬ 
tives  of  second-order  coefficients  omitting  dependence  of  elastic 
coefficients  on  shear  strains  [23]: 

CafySa<l>  =  ~  {BI'i)\{dCapxi/dp)Seti>+  {dCE,j>nfj dp) Sxg 

+  U>CxSerf/dp)Srlfi]  (12) 

Dislocations  are  addressed  in  Table  4:  Full  1/2(110)  and  partial 
1/4(110)  dislocations  are  considered.  Line  energy  W  (full  and 
partial  dislocations)  and  Peierls-Nabarro  stress  tpn  (full  disloca¬ 
tions  only)  are  calculated  analytically  [16]  via 


IV  =  \Kb2l(4TT)]log(R/R0)  =  Kb2,  rPN  =  2 K  exp[-  2t t(KI/1) 

X(d/b )]  (13) 

where  R  is  the  distance  from  the  dislocation  line,  R0  is  the  core 
radius,  (l  is  the  appropriate  shear  elastic  constant  (see  Appendix), 

d  is  the  spacing  between  slip  planes,  and  k  depends  on  the 
second-order  elastic  coefficients  and  the  geometry  of  the  disloca¬ 
tion  line.  Partial  dislocations  on  {111}  planes  have  been  observed 
in  dynamically  compressed  AlON  [4],  Slip  on  {110}  planes  is  also 
considered  in  Table  4  since  near  room  temperatures,  Veyssiere  et 
al.  [8]  found  slip  resistance  on  {110}  planes  to  be  lower  than  that 
on  {111}  in  Mg0-nAl203  spinel  (n~  1.1)  single  crystals.  On  the 
other  hand,  Mg0-nAl203  spinel  (n~  3.5)  exhibited  a  lower  re- 


Table  3  Elastic  properties 


Parameter 

Value 

Definition 

Reference 

Cu 

377  GPa 

Second-order  elastic  constants 

[10] 

C\ 2 

133  GPa 

C44 

125  GPa 

Ao 

124  GPa 

Shear  modulus  (l/5)(C11-C12)  +  (3/5)C44 

Voigt  average 

Bo 

214  GPa 

Isothermal  bulk  modulus  (1 /3)(C11  +  2C12) 

Bv 0 

216  GPa 

Isentropic  bulk  modulus 

Eq.  (6) 

dfil  dp 

0.95 

Pressure  derivative  of  shear  modulus 

[13] 

dBIdp 

4.2 

Pressure  derivative  of  bulk  modulus 

[13] 

dpi  38 

-0.012  GPa/K 

Temperature  derivative  of  shear  modulus 

[13] 

3BI38 

-0.017  GPa/K 

Temperature  derivative  of  bulk  modulus 

[13] 

Table  4  Theoretical  dislocation  properties 


Plane 

Burgers 

b 

(A) 

d 

(A) 

Line  and  type3 

k 

(GPa) 

A  . 
(GPa)b 

W 

(eV/A)c 

rPN 

(GPa) 

{in} 

1/2[1 10] 

5.62 

4.59 

[110]  screw 

123 

124 

24.4 

1.49 

[101]  mixed 

156 

30.1 

0.49 

Isotropic  edge 

167 

32.9 

0.34 

1/4[1 10] 

2.81 

4.59 

[110]  screw 

123 

6.1 

- 

[101]  mixed 

156 

7.7 

- 

Isotropic  edge 

167 

8.2 

{110} 

1/2[1 10] 

5.62 

5.62 

[110]  screw 

123 

122 

24.4 

0.43 

Isotropic  edge 

167 

32.9 

6.2  X  10-2 

1/4[1 10] 

2.81 

5.62 

[110]  screw 

123 

6.1 

- 

Isotropic  edge 

167 

8.2 

- 

{100} 

1/2[1 10] 

5.62 

7.95 

[110]  screw 

123 

125 

24.4 

3.8  x  10-2 

Isotropic  edge 

167 

32.9 

2.4  X  10-3 

1 /4[ 1 1 0] 

2.81 

7.95 

[110]  screw 

123 

6.1 

- 

Isotropic  edge 

167 

8.2 

- 

dAnisotropic  solutions  for  screw  and  mixed  only;  isotropic  solutions  for  edge  valid  for  any  line  direction. 
bFormulae  listed  in  Appendix. 

0 Assuming  log(/?//?0)=4/77\ 
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Table  5  Theoretical  shear  fracture  and  twinning  properties 


Plane 

Burgers 

b 

(A) 

d 

(A) 

A 

(GPa) 

tF 

(GPa) 

y 

(J  /  m2) 

tt 

(GPa) 

{in} 

1/2[1 10] 

5.62 

4.59 

124 

24.2 

_ 

_ 

1/12[1 12] 

1.62 

- 

- 

- 

0.104“ 

0.64 

1/12[112] 

1.62 

- 

- 

- 

0.140b 

0.86 

{110} 

1/2[1 10] 

5.62 

5.62 

122 

19.4 

- 

- 

{100} 

1/2[1 10] 

5.62 

7.95 

125 

14.1 

" 

- 

‘'Screw  dislocation. 
bEdge  dislocation. 


sistance  to  slip  on  {1 1 1}  planes  than  on  {110}  or  {100}  [8].  Experi¬ 
ments  of  Veyssiere  et  al.  [8]  involved  the  compression  of  single 
crystals  under  confining  pressure.  In  contrast,  studies  of  Knoop 
(indentation)  microhardness  profiles  of  MgO-Al203  and 
Mg0-3A1203  spinel  single  crystals  [36]  suggest  that  (1 10){1 11} 
are  active  slip  systems  at  room  temperature  although  such  experi¬ 
ments  suggest  secondary  slip  on  (110){110}  could  be  possible  in 
the  latter  (nonstoichiometric)  spinel.  In  recovered  samples  of 
shock-loaded  polycrystalline  MgO-Al203,  full  1/2(110){111} 
dislocations  and  {111}  twin  lamellae  have  been  observed  via  trans¬ 
mission  electron  microscopy  [37],  For  completeness,  slip  on  {100} 
is  also  addressed  in  Table  4  since  (110){100}  slip  has  been  ob¬ 
served  in  spinel  crystals,  albeit  with  higher  room  temperature  re¬ 
sistance  than  ( 1 10){  1 10}  slip  [8],  From  Table  4,  analytical  Peierls- 
Nabarro  stresses  are  greater  for  slip  on  {111}  planes  than  for  slip 
on  {110}  planes,  in  qualitative  agreement  with  trends  reported  for 
Mg0-nAl203  spinel  (n~  1.1)  [8],  However,  the  Peierls-Nabarro 
stress  for  slip  on  {100}  is  much  lower  than  that  for  slip  on  {111},  in 
disagreement  with  experimental  shear  strengths  reported  for  other 
spinels  [8,36].  Dislocation  core  energies  often  obtained  via  atomic 
calculations  are  presently  unknown  for  AlON  and  are  not  ad¬ 
dressed  in  Table  4.  Properties  pertaining  to  twinning  and  shear 
fracture  are  listed  in  Table  5.  Twinning  stress  and  theoretical  shear 
fracture  strength  are  estimated,  respectively,  as  [2,17,23] 

rT=ylb ,  tf  =  jlblilird)  (14) 

where  b  in  the  first  of  Eq.  (14)  is  the  magnitude  of  the  Burgers 
vector  of  the  partial  dislocation  associated  with  twinning  and  b  in 
the  second  of  Eq.  (14)  is  the  magnitude  of  a  full  Burgers  vector  on 
the  cleavage  plane.  Twinning  is  considered  to  occur  on  {111} 
planes  with  corresponding  partial  Burgers  vector  1/12(112)  [38], 
Stacking  fault  energy  y  is  estimated  with  an  isotropic  linear  elastic 
approximation  [9]: 

y=  /ib2/(2nReq)[cos2  (f>+  sin2  </>/(l  -  n)]  (15) 

where  4>  is  the  angle  between  the  Burgers  vector  and  dislocation 
line,  the  equilibrium  spacing  Req  between  partial  dislocations  is  15 
nm  [4],  and  v=(3B-2p)/(6B+2p)  is  the  Poisson’s  ratio.  Com¬ 
puted  {111}  stacking  fault  energy  of  AlON  in  Table  5  is  compa¬ 
rable  to  that  observed  for  MgAl204  spinel  [9].  Properties  in  Tables 
4  and  5  serve  as  guidance  for  continuum  model  parameters  per¬ 
taining  to  defects  entering  Eqs.  (4)  and  (11),  as  discussed  in  Sec. 
4.  Peierls-Nabarro  stresses  have  been  used  elsewhere  to  explain 
plastic  slip  resistance  in  ceramic  single  crystals  [1,2]  and  theoret¬ 
ical  strengths  have  been  considered  elsewhere  for  use  as  failure 
criteria  for  other  ceramics  subjected  to  impact  loading  [2,20]. 

4  Model  Application:  High-Pressure  Loading  of  Alu¬ 
minum  Oxynitride 

The  model  outlined  in  Sec.  2  is  applied  to  AlON  [39]single-  and 
polycrystals  deformed  under  high-pressure  loading  as  occurring  in 
plate  impact  experiments  [12,24-27].  Because  data  for  the  impact 
response  of  individual  single  crystals  are  not  available,  polycrys¬ 
talline  data  are  used  to  guide  selection  of  certain  model  param- 

011013-4  /  Vol.  78,  JANUARY  2011 


eters  for  slip,  twinning,  pore  collapse,  and  shear  fracture.  For  in¬ 
stances  wherein  sufficient  data  are  not  available  to  enable  unique 
parameter  selection  (e.g.,  initial  resistances  for  slip  and  twinning 
on  various  families  of  crystallographic  planes),  several  parameter 
sets  are  explored. 

4.1  Slip  and  Twinning  Parameters.  The  following  relation¬ 
ships  are  used  for  slip  and  twin  resistances  entering  flow  rules  in 
Eq.  (11): 

7c=(7co  +  A'7c)(mW>  tc=(tco  +  ^tc)(aW  (16) 

&Tic=h(^L  y+S/vj.  Af,c=/t|2  y+2/yj  07) 

In  Eq.  (16),  t}0  and  t}0  are  the  constant  initial  values  of  slip 
and  twin  resistances  that  may  differ  among  families  of  slip  and 
twin  systems,  At},  and  At},  account  for  the  changes  in  slip  and 
twinning  resistance  with  the  changes  in  defect  densities  at  stan¬ 
dard  temperature  and  pressure,  fj,0  is  the  shear  modulus  at  stan¬ 
dard  temperature  and  pressure  (the  constant  Voigt  average  in  Table 
3),  and  p  is  the  effective  shear  modulus  that  depends,  like  bulk 
modulus  B,  linearly  on  temperature  and  pressure 

fi  =  p0  +  p(dpldp)  +  (8-  80)(dpldd),  B  =  B0  +  p(8B/dp)  +  (8 
-  80)(dB/dO)  (18) 

In  the  first  of  Eq.  (16),  scaling  of  the  slip  resistance  by  p!  pQ 
follows  the  assertion  that  slip  resistance  should  be  proportional  to 
the  instantaneous  shear  modulus  [40,41]  and  results  in  relatively 
mild  influences  of  pressure  and  temperature  on  yield  strength.  The 
variation  of  yield  stress  with  temperature  could  conceivably  be 
more  substantial,  e.g.,  represented  by  Arrhenius-type  kinetics  or 
temperature-dependent  dislocation  drag,  but  sufficient  data  are  not 
available  for  AlON  to  justify  a  more  complicated  relation.  Twin¬ 
ning  resistance  follows  an  analogous  prescription  in  the  second  of 
Eq.  (16).  In  Eq.  (17),  h  is  a  constant  hardening  modulus  and  all 
slip  and  twin  systems  harden  at  the  same  rate,  the  simplest  real¬ 
istic  assumption  possible  (i.e.,  equal  self-  and  latent-hardening). 
The  model  in  Eq.  (17)  represents  the  following  general  physics 
[17,42]:  Accumulated  dislocations  and  twin  boundaries  act  as  bar¬ 
riers  to  slip  and  further  twinning.  The  relationship  between  accu¬ 
mulated  dislocation  density  and  change  in  slip  resistance  is  as¬ 
sumed  to  follow  the  usual  Taylor-type  law  [2,43-45]  deemed 
applicable  for  spinel  ceramics  [9]: 


where  a  for  most  materials  lies  between  0.1  and  1.0  [9].  The  area 
per  unit  volume  of  twin  boundaries  obtained  from  simple  geomet¬ 
ric  considerations  is 

777-  =  2^/-,/tj.=  2/r/;j.  (20) 

i 

with  tT  as  the  thickness  of  twin  lamellae  characteristic  of  twin 
structures  observed  in  AlON  deformed  at  high  rates  [4,7], 

Parameters  entering  the  single  crystal  constitutive  model  are 
listed  in  Table  6.  Hardening  modulus  h  is  chosen  to  simulta- 
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Table  6  Parameters  for  a  single  crystal  plasticity  model 


Parameter 

Value 

Definition 

h 

1.0  GPa 

Hardening  modulus 

To 

1.0/s 

Reference  strain  rate 

m  and  n 

10 

Strain  rate  sensitivity 

y 

V2/2 

Twinning  shear 

b 

5.62  A 

Burgers  magnitude 

k 

123-0.012(0- 0O)  GPa 

Dislocation  energy  factor 

y 

0.14  J/m2 

Stacking  fault  energy 

a 

0.5 

Proportionality  constant 

tT 

10  /im 

Characteristic  twin  thickness 

neously  produce  an  accurate  representation  of  the  shear  stress  re¬ 
sponse  of  polycrystalline  AlON  subjected  to  uniaxial  strain  load¬ 
ing  and  a  physically  reasonable  dislocation  density  in  Eq.  (19),  as 
discussed  later  in  Sec.  4.2.  Rate  sensitivity  parameters  follow 
from  experimental  observations  on  alumina  (a,-Al203)  [46],  Note 
that  AlON,  stoichiometric  ally,  is  related  to  both  alumina  and  alu¬ 
minum  nitride  [7,11],  so  use  of  alumina  data  to  justify  selection  of 
m  may  be  a  reasonable  first  assumption  in  the  absence  of  test  data 
for  AlON  or  spinel.  In  many  crystalline  solids,  inverse  rate  sensi¬ 
tivity  1/m  increases  with  loading  rate;  the  value  chosen  here  for 
impact  loading  (l/m  =  0. 1)  may  not  apply  at  lower  strain  rates. 
Coefficient  yQ  provides  dimensional  consistency.  Twinning  shear 
y  is  standard  for  a  compound  twin  in  spinel  [38,47],  Representa¬ 
tively,  b  is  taken  as  the  full  Burgers  vector  for  (110)  slip,  K  is 
chosen  as  the  anisotropic  energy  factor  for  a  screw  dislocation 
(Table  4)  with  a  temperature  dependence  equal  to  that  of  the  shear 
modulus  (Table  3),  stacking  fault  energy  y  is  that  for  a  twinning 
edge  partial  (Table  5),  a  in  Eq.  (19)  is  assigned  a  typical  value  of 
0.5,  and  the  characteristic  twin  thickness  is  10  pm  [7]  on  the 
order  of  1/10  of  a  typical  grain  diameter. 

The  following  slip  systems  are  incorporated  in  subsequent  cal¬ 
culations:  {1 1 1}(  1 10),  {1 10}(1 10),  and  {100}(110),  all  which  have 
been  observed  in  spinel  ceramics  [8,9].  Twinning  systems  are  of 
type  {111}(112)  [38],  Four  sets  of  values  for  initial  hardness  of 
slip/twin  systems  in  Eq.  (16)  are  explored,  as  listed  in  Table  7. 
The  first  of  set  assumes  that  all  the  slip  and  twinning  resistances 
are  equal,  the  simplest  assumption  incorporating  all  possible 
mechanisms.  The  second  set  assumes  that  the  only  significant  in¬ 
elastic  deformation  mode  is  { 1 1 1  }<  1 1 0)  slip,  in  agreement  with 
indentation  anisotropy  profiles  for  Mg0-Al203  spinel  single 
crystals  [36],  which  do  not  seem  to  twin  readily  [9]  although  [111] 
twin  lamellae  have  been  observed  in  shock  recovered  specimens 
[37],  The  third  set  assumes  that  a  ratio  of  initial  resistances  for 
shear  on  (111}:{11 0}:{  1 00}  planes  of  2:1:2,  motivated  by  experi¬ 
mental  observations  for  Mg0-uAl203  spinel  (n~  1.1)  [8]  at  or 
below  400  °C.  The  fourth  set  assumes  that  the  ratio  of  initial  twin¬ 
ning  resistance  to  initial  slip  resistance  on  all  systems  is  1:2.  Ini¬ 
tial  resistances  for  slip  are  comparable  to  Peierls-Nabarro  stresses 
for  screw,  mixed,  or  edge  dislocations  on  {111}  and  {110}  planes 
computed  from  Eq.  (13)  and  listed  in  the  rightmost  column  of 


Table  4;  initial  resistances  for  twinning  are  comparable  to  theoret¬ 
ical  estimates  computed  from  Eq.  (14)  and  listed  in  the  rightmost 
column  of  Table  5. 

Twinning  induces  a  1 80  deg  reflection  or  rotation  of  the  lattice 
vectors  about  the  habit  plane  normal  or  about  the  twin  axis 
[42,47],  For  compound  twins  in  centrosymmetric  crystals,  such  an 
orthogonal  transformation  for  twin  system  j  can  be  represented  by 
matrix  Qlap=:2mlanTlp- Sap,  a  180  deg  rotation  about  habit  plane 
unit  normal  mJa.  Re-orientation  of  the  elastic  modulus  tensor  due 
to  twinning  [2]  is  not  required  in  the  present  study  because  of  the 
low  elastic  anisotropy  of  AlON  (Zener’s  ratio  2C44/(C’I1-Ci2) 
=  1.02—1)  and  lattice  orientation  does  not  affect  the  thermal  ex¬ 
pansion  coefficient  for  a  cubic  crystal.  Twinning  is  unidirectional 
in  Eq.  (11);  when  the  resolved  shear  stress  is  negative  on  the  habit 
plane  in  the  direction  of  positive  shear,  inelastic  deformation  on 
{111}  planes  can  result  only  from  glide  of  partial  dislocations.  In  a 
previous  model  [2],  the  summed  contribution  from  slip  rates  in  the 
first  of  Eq.  (3)  was  weighted  by  factor  1  —fT  and  new  secondary 
slip  systems  could  become  active  within  the  twins.  In  the  present 
model,  factor  1  —fT  is  omitted  but  an  offsetting  resistance  to  slip  is 
reflected  by  hardening  due  to  twinning  in  Eq.  (17).  For  simplicity 
and  because  of  insufficient  microscopy  data  for  AlON  to  justify  a 
more  complicated  theory,  new  slip  systems  are  not  activated 
within  twins  in  the  present  application  and  successive  twins  of  a 
different  orientation  within  a  given  twin  (including  reversal  or 
detwinning)  are  not  modeled  explicitly.  In  calculations  that  follow, 
the  volume  fraction  of  twins  in  any  crystal  always  remains  less 
than  20%  (preferred  {111}  twinning  in  Table  7)  and  usually  re¬ 
mains  less  than  5%  (all  other  cases  in  Table  7).  This  model  might 
be  inappropriate  for  conditions  wherein  the  total  twin  volume 
fraction  approaches  unity. 

4.2  Polycrystal  Modeling.  The  porosity  of  a  polycrystal  can 
significantly  affect  its  mechanical  response.  To  first  address  the 
effect  of  porosity  independently  from  slip  or  twinning,  the  re¬ 
sponse  of  an  aggregate  of  grains  subject  to  hydrostatic  pressure  is 
studied.  For  a  cubic  crystal,  the  deviatoric  and  volumetric  stress- 
strain  responses  are  decoupled  and  the  stress  in  a  cubic  crystal 
subject  to  spherical  compression  is  hydrostatic  pressure,  indepen¬ 
dent  of  the  crystal’s  lattice  orientation.  For  a  prescribed  change  in 
volume,  the  average  pressure  p  supported  by  a  polycrystal  de¬ 
pends  on  the  aggregate’s  effective  bulk  modulus,  porosity,  and 
temperature.  The  influence  of  porosity  on  effective  bulk  modulus 
B  and  effective  shear  modulus  p  of  the  polycrystalline  aggregate 
is  represented  here  using  the  self-consistent  model  of  Mackenzie 
[48] 

B  =  4pB{l  -  d)l[4p  +  3 ##],  p  =  p\\  -  5tf(3 B  +  4p)H9 B 
+  8/r)]  (21) 

where  #  is  the  porosity  and  p  and  B  are  the  effective  shear  and 
bulk  moduli  of  the  fully  dense  crystal  that  depend  on  pressure  and 
temperature,  as  in  Eq.  (18).  Above  a  threshold  pressure,  pores  may 
collapse  irreversibly.  This  phenomenon,  common  in  brittle  ceram¬ 
ics,  minerals,  and  geologic  materials  [18,22],  is  modeled  via  [19] 


Table  7  Initial  resistances  for  slip  and  twin  systems  in  a  crystal  plasticity  model 


Slip  or  twin 
system 

Tco 

all  systems 
equal 
(GPa) 

TCO 

exclusive 
{111}  slip 
(GPa) 

TC0 

preferred 
{110}  slip 
(GPa) 

TC0 

preferred 
{111}  twinning 
(GPa) 

{11 1}<1 10)  slip 

1.00 

0.75 

1.10 

1.50 

{110}(110>  slip 

1.00 

00 

0.55 

1.50 

{look  110)  slip 

1.00 

OO 

1.10 

1.50 

{111}(112>  twin 

1.00 

OO 

1.10 

0.75 
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Fig.  1  Hydrostatic  pressure  for  polycrystalline  AION:  experi¬ 
ments  [25]  and  model  predictions 


sj}jRSp[mp-p0)H(Sp)]  [*>0] 

[o  [tf=0] 

where  S-&  is  an  incremental  reduction  in  porosity  commensurate 
with  average  pressure  increase  Sp  supported  by  the  aggregate,  R  is 
a  negative  constant  dictating  the  rate  of  pore  collapse,  H(.)  is  a 
Heaviside’s  step  function,  and  p0  is  a  threshold  pressure. 

Model  results  are  shown  in  Fig.  1  for  AION  subjected  to  spheri¬ 
cal  compression.  Experimental  data  are  from  impact  experiments 
[25]  reproduced  in  Ref.  [12].  In  the  application  of  the  model  cor¬ 
responding  to  Fig.  1,  adiabatic  conditions  are  assumed  and  ther¬ 
moelastic  heating  is  taken  into  account  (fifth  term  on  the  right  in 
Eq.  (9))  but  no  slip  or  twinning  occurs  because  no  shear  stresses 
are  generated  during  spherical  compression.  The  average  pressure 
is  p  and  the  total  volume  change  of  the  aggregate  (crystals  and 
pores)  is  measured  by  Jacobian 


crystals  follow  the  Taylor  assumption  [44,49,50]:  All  grains  in  the 
aggregate  are  subject  to  same  deformation  gradient  components 
Fa  A-  This  assumption  is  deemed  more  appropriate  than  other 
simple  averaging  schemes  (e.g.,  those  that  assume  constant  stress 
in  all  grains)  for  highly  confined  loading  conditions  [50],  For 
uniaxial  strain,  the  deformation  gradient  in  matrix  form  is 


1  -et  0  0 


[F]  = 


0 

0 


1  0 
0  1 


(24) 


where  e  is  the  applied  strain  rate  positive  for  compression.  The 
rate  used  in  all  calculations  that  follow  is  e=105/s,  representative 
of  plate  impact  [21].  For  the  present  set  of  boundary  conditions, 
elastic  volume  change  JE  is  identical  in  each  grain  and  Eq.  (23) 
again  describes  the  total  volume  change  of  the  porous  aggregate. 
Deviatoric  elastic  and  plastic  deformations  differ  among  grains 
with  different  initial  orientations.  Average  and  local  deviatoric 
Cauchy  stress  components  are  defined,  respectively,  as 


® ab  G 'ab  PBab 


(25) 


with  pi  computed  from  Eq.  (21)  using  the  porosity,  the  average 
temperature,  and  the  average  pressure  in  the  grains  to  determine 
elastic  coefficients  of  the  matrix  in  Eq.  (18).  Average  temperature, 
twin  volume  fraction,  and  dislocation  density  in  the  grains  are 


0=(l/AOl>,  /•/  =  (] //V)2/r:  and  pD  =  (l/N)^pD 

N  N  N 

(26) 

The  average  pressure  for  the  porous  polycrystalline  aggregate  is 


y  =  ./£[i -(#„-#)]  (23) 

where  -&0  is  the  initial  porosity  and  JE  is  the  elastic  compression 
that  is  identical  in  each  grain.  Four  model  curves  are  shown  with 
the  test  data.  The  most  accurate  fit  is  provided  by  the  curve  la¬ 
beled  “pore  crush,”  which  uses  Eqs.  (21)  and  (22),  with  param¬ 
eters  for  Eq.  (22)  listed  in  Table  8.  Initial  porosity  -&0  is  consistent 
with  the  initial  mass  density  of  specimens  used  in  experiments 
[25].  The  curve  labeled  “fully  dense”  corresponds  to  predictions 
for  a  crystal  with  null  porosity.  The  curve  labeled  “fixed  porosity” 
assumes  no  pore  collapse  (R  =  0  in  Eq.  (22))  but  allows  for  a 
reduction  in  bulk  modulus  via  a  fixed  value  of  rl0  in  Eq.  (21).  Also 
shown  is  a  curve  for  uniaxial  strain  loading  conditions  that  will  be 
discussed  later;  this  “uniaxial  strain”  curve  is  obtained  using  Eqs. 
(21)  and  (22)  and  is  nearly  identical  to  the  response  for  spherical 
compression.  The  regime,  where  the  slope  of  pressure  versus  vol¬ 
ume  is  reduced,  corresponds  to  contributions  from  collapsing 
pores.  From  Fig.  1,  when  pore  collapse  is  omitted,  average  pres¬ 
sure  is  severely  overestimated  for  compression  greater  than  4% 
(i.e.,  for  7<0.96). 

Considered  next  is  the  response  of  an  aggregate  of  N  randomly 
oriented  crystals  of  equal  volume  subjected  to  uniaxial  strain.  The 
aggregate  also  contains  initial  porosity  U0,  which  can  be  concen¬ 
trated  at  boundaries  between  grains.  Results  that  follow  for  poly- 


Table  8  Parameters  for  polycrystal  pore  collapse  and  shear 
fracture  models 


Parameter 

Value 

Definition 

0.033 

Initial  porosity 

Po 

10.0  GPa 

Initiation  pressure  for  pore  collapse 

R 

-0.005 /GPa 

Pore  collapse  coefficient 

3.5  GPa 

Macroscopic  shear  fracture  strength 

X 

3.25  X10“5  GPa-s 

Post-fracture  dynamic  viscosity 

p  = - (27) 

NBn 

where  B  is  computed  from  Eq.  (21).  The  maximum  shear  stress 
for  the  aggregate  is  defined  as 

T=(^Lx~BmJ/2  (28) 

where  CT^ax  and  amin  are  the  maximum  and  minimum  principle 
values  of  the  first  of  Eq.  (25).  Under  confined  uniaxial  strain  com¬ 
pression  Eq.  (24),  the  present  material  (consisting  of  cubic  crys¬ 
tals  and  pores)  experiences  no  tensile  stress.  However,  above 
some  threshold  maximum  shear  stress  tf,  the  material  may  fail 
due  to  cleavage  or  grain  boundary  fracture.  After  such  failure 
occurs,  the  material  consists  of  numerous  fragmented  grains  and  is 
represented  as  an  isotropic,  compressible,  and  viscous  medium. 
Considered  in  some  calculations  that  follow  is  a  shear  failure 
model  that  takes  effect  for  providing  the  following  rela¬ 

tions  for  average  deviatoric  stress  and  temperature  rise  in  the  post¬ 
fracture  regime: 

a'ab  =  2  XD'ab,  d=JEa'abD'j(p0cv)  (29) 

The  deviatoric  deformation  rate  is  D'ab=(Lab+Lba)l2 
-SabLccl?>.  The  post-fracture  dynamic  viscosity  is  \  and  is  as¬ 
sumed  constant  for  simplicity  in  Table  8  although  a  more  general 
treatment  would  permit  possible  dependence  of  sliding  friction  of 
fragmented  grains  on  pressure  and,  perhaps,  temperature  and  load¬ 
ing  rate  [22].  For  rS  tf,  pressure  is  computed  in  same  way  as  in 
Eq.  (27).  The  second  of  Eq.  (29)  represents  the  temperature  rise 
from  dissipation  due  to  shear  friction  of  fractured  grains  sliding 
past  one  another.  Numerical  algorithms  for  the  present 
thermoelastic-plastic  calculations  extend  those  used  in  previous 
work  for  large,  high-rate  deformations  of  cubic  crystals  [31,44,45] 
to  account  for  aforementioned  slip,  twinning,  pore  collapse,  and 
shear  fracture  mechanisms  thought  relevant  for  AION 
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Fig.  2  Axial  stress  for  polycrystalline  AION  in  uniaxial  strain:  (a)  experiments  [24-27]  and  model 
predictions  (equal  slip  and  twin  resistances)  and  ( b )  model  predictions  for  various  slip/twin/ 
fracture  parameters  and  elasticity  models  (upper  four  curves  in  (b)  are  nearly  indistinguishable) 


polycrystals. 

Stress-strain  data  and  model  predictions  are  shown  in  Fig.  2. 
Axial  stress  supported  by  the  aggregate  is 

^=-(°h-p)  (30) 

positive  in  compression.  Predictions  shown  in  Fig.  2(a)  assume 
that  all  slip/twin  systems  have  equal  hardness  (column  2  of  Table 
7).  Predictions  in  Fig.  2(b)  account  for  the  four  sets  of  initial 
slip/twin  system  hardness  values  listed  in  columns  of  Table  7  and 
for  the  shear  fracture  model  of  Eq.  (29).  Predictions  for  the  latter 
(shear  fracture)  model  also  assume  equal  hardness  among  slip/ 
twin  systems  for  t<tf.  Also  shown  in  Fig.  2(b)  are  results  la¬ 
beled  “linear  elasticity,”  which  are  obtained  using  the  same  model 
parameters  as  “all  systems  equal”  with  the  exception  that 
pressure-  and  temperature-derivatives  of  elastic  moduli  are  set  to 
zero.  All  model  results  that  incorporate  nonlinear  elasticity  are 
nearly  indistinguishable  in  Fig.  2(b)  and  thus  provide  an  accurate 
fit  to  the  plate  impact  data  [24-27]  in  Fig.  2(a).  Model  results 
incorporating  linear  elasticity  noticeably  underpredict  the  axial 
stress  for  large  compressions,  i.e.,  for  7<0.97.  Notice  that  pores 
are  crushed  according  to  Eq.  (22);  their  contribution  to  volume 
change,  deviatoric  stress,  and  pressure  is  incorporated  in  Eqs. 
(23),  (25),  and  (27),  respectively.  It  was  determined  that  N=  100 
(i.e.,  100  random  grain  orientations)  was  sufficient  to  ensure  a 
macroscopically  isotropic  response  insensitive  to  further  increases 
in  N. 

Shear  stresses  supported  by  the  aggregate  are  shown  in  Fig.  3; 
predictions  are  computed  via  Eq.  (28).  All  predictions  in  Fig.  3 
provide  a  reasonable  description  of  the  experimental  data  [12,25], 
which  exhibit  notable  scatter.  As  summarized  elsewhere  [12],  ex¬ 


J 

Fig.  3  Maximum  shear  stress  for  polycrystalline  AION  in 
uniaxial  strain:  experiments  [25]  and  predictions  for  slip/twin/ 
fracture  model  parameters 

Journal  of  Applied  Mechanics 


periments  by  other  researchers  have  indicated  that  the  shear 
strength  of  AION  polycrystals  shocked  above  the  Hugoniot  elastic 
limit  (HEL)  typically  falls  in  the  range  3. 3-3.9  GPa,  in  general 
agreement  with  model  predictions  in  Fig.  3.  The  experimental 
HEL  corresponds  to  a~  10-12  GPa  or  7  —  0.97  [12,26,27],  Ini¬ 
tial  resistances  to  slip  or  twinning  used  in  the  model  (Table  7)  are 
~1  GPa,  of  the  order  of  Peierls-Nabarro  stresses  for  {111}  and 
{110}  slip  of  screw  dislocations  (Table  4)  and  comparable  to  shear 
strengths  for  slip  in  other  spinel  crystals  (—1-2  GPa)  [8]  but  are 
an  order  of  magnitude  smaller  than  theoretical  strengths  predicted 
for  AION  (Table  5).  At  larger  compressive  strains,  predicted  shear 
stresses  will  continue  to  increase  for  all  cases  except  the  model 
incorporating  shear  fracture  for  which  shear  stress  increases  little 
for  7<0.93.  Plate  impact  experiments  on  AION  have  demon¬ 
strated  a  reduction  in  spall  strength  from  1.7  GPa  to  0.14  GPa  as 
the  impact  stress  is  increased  from  4.8  GPa  to  9.5  GPa  [24,27], 
indicating  that  some  local  shear-induced  damage  or  fracture  pro¬ 
cesses  may  take  place  at  impact  pressures  below  9.5  GPa.  It  is 
thus  likely  that  a  model  incorporating  shear  fracture  (e.g.,  via  Eq. 
(29))  is  physically  realistic  for  AION  polycrystals  subject  to  im¬ 
pact  stresses  on  the  order  of  the  HEL.  Hardening  law  for  slip  Eq. 
(17)  could  be  modified  so  that  strain  hardening  ceases  above  some 
saturation  stress  [29,44,45];  however,  from  the  preceding  argu¬ 
ments,  it  may  be  reasonable  to  assume  that  during  monotonic 
impact  loading,  hardening  and  accumulation  of  dislocations  by 
Eq.  (19)  and  twin  boundaries  by  Eq.  (20)  proceed  until  fracture. 
Intragranular  fracture  may  be  initiated  as  a  result  of  local  stress 
concentrations  due  to  interactions  of  dislocations  and  twin  bound¬ 
aries  [42]. 

Predictions  for  average  temperature  rise  during  adiabatic 
uniaxial  straining  of  the  aggregate  are  shown  in  Fig.  4.  wherein 
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Fig.  4  Average  temperature  for  poiycrystalline  AION  in 
uniaxial  strain:  model  predictions 
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Fig.  5  Predicted  defect  densities  for  polycrystalline  AION  in  uniaxial  strain:  (a)  total  twin  volume 
fraction  and  ( b )  dislocation  density 


the  local  temperature  of  each  grain  is  updated  according  to  Eq.  (9) 
with  r=  0  and  V2/?=0.  The  different  initial  slip/twin  hardness  val¬ 
ues  and  the  use  of  shear  fracture  model  all  result  in  similar  aver¬ 
age  temperature  predictions.  In  hydrostatic  compression  calcula¬ 
tions  corresponding  to  Fig.  1,  the  temperature  rise  resulting  from 
thermoelasticity  alone  (fifth  term  on  the  right  of  Eq.  (9))  is 
~20  K  at  .7=0.87.  For  uniaxial  strain  corresponding  to  Fig.  4  at 
7=0.87,  the  adiabatic  temperature  rise  resulting  from  plastic  dis¬ 
sipation  less  the  energy  storage  associated  with  lattice  defects  as¬ 
sociated  with  the  rightmost  two  terms  in  Eq.  (9)  is  ~100  K.  Pos¬ 
sible  temperature  rise  in  the  material  associated  explicitly  with 
pore  collapse  (e.g.,  adiabatic  heating  in  highly  plastically  de¬ 
formed  regions  surrounding  voids  and  heat  transfer  between  crys¬ 
tals  and  compressed  gas/vapor  within  voids)  is  omitted  in  the 
present  analysis,  as  is  that  resulting  from  any  additional  local 
dissipation  associated  with  passage  of  a  shock  front. 

Predicted  twin  volume  fractions  are  shown  in  Fig.  5(n).  When 
the  critical  shear  stress  required  for  twinning  is  reduced,  as  dem¬ 


onstrated  in  results  for  “preferred  {111}  twinning,”  the  volume 
fraction  of  twins  produced  during  the  deformation  process  in¬ 
creases.  When  {110}  slip  is  favored,  on  the  other  hand,  the  fraction 
of  twins  produced  decreases  relative  to  results  in  which  all  slip 
and  twin  systems  have  equal  strength.  During  all  calculations,  the 
maximum  volume  fraction  of  twins  in  any  grain  always  remained 
less  than  0.2.  For  calculations  incorporating  exclusive  {111}  slip, 
twinning  is  prohibited  so  the  volume  fraction  of  twins  for  that 
case  is  zero.  The  present  results  suggest  that  experimental  mea¬ 
surements  of  twin  volume  fractions  in  recovered  samples  should 
provide  insight  into  possible  preference  of  slip  or  twinning  in 
AION  deformed  at  high  rates  under  confining  pressure.  Figure 
5(b)  shows  predicted  dislocation  densities.  Magnitudes  are  com¬ 
parable  to  those  observed  for  other  spinel  crystals  undergoing 
plastic  deformation  [9].  Experimental  and  theoretical  studies  have 
suggested  that  local  residual  stress  fields  of  dislocations  may  im¬ 
part  a  (small)  volume  change  in  crystals,  which  can  be  estimated 
as  [2,31,32] 


SV/1-- 


■Iv2 


(1  ~v)B 

1  /  dfJL  yU. 

yU.  \  dp  B 


(H 


:  -  2v+  2 it  i  dpi  fjL 
(1  -  v)fjL  \dp  B  / 


W 


W  (edge  dislocations) 
(screw  dislocations) 


(31) 


where  SV/l  is  the  volume  change  per  unit  length  of  dislocation 
line  and  VP  is  the  energy  per  unit  length  (Table  4).  From  Eq.  (31) 
and  Tables  3  and  4,  SVIl= 3.3  X  10“19  m2  for  1/2(110)  edge  dis¬ 
locations  and  SV/!=  1.2X10-19  m2  for  1/2(110)  screw  disloca¬ 
tions,  predictions  comparable  to  those  for  other  ceramics  [2,23], 
For  a  dislocation  density  of  pfl=1013  m-2,  the  volumetric  strain 
associated  with  such  (edge)  dislocations  is  pDX  SVI /=3.3X  10-6, 
negligible  in  comparison  to  applied  volumetric  strains  on  the  or¬ 
der  of  0.01.  Thus,  omission  of  volume  changes  associated  with 
dislocations,  as  indicated  in  Eqs.  (1)— (3),  is  a  reasonable  modeling 
assumption  for  the  present  loading  conditions. 

4.3  Single  Crystal  Predictions.  Investigated  next  are  re¬ 
sponses  of  single  crystals  of  several  orientations  subjected  to 
uniaxial  strain  Eq.  (24)  at  a  fixed  rate  of  e=105/s  and  again  as¬ 
suming  adiabatic  conditions.  Crystals  are  treated  in  Sec.  4.3  as 
fully  dense  (i.e.,  null  porosity)  and  failure  associated  with  fracture 
at  grain  boundaries  addressed  via  Eq.  (29)  for  polycrystals  is  not 
considered  for  single  crystals.  Two  initial  lattice  orientations  are 
considered.  The  first,  labeled  [100],  corresponds  to  the  alignment 
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of  the  crystal’s  cube  axes  with  the  global  coordinate  system,  so 
that  compression  in  Eq.  (24)  is  normal  to  the  (100)  plane.  The 
second  orientation,  labeled  [110],  corresponds  to  Goss  orientation 
of  texture  analysis  [51]  with  Bunge’s  Euler  angles  (<Pi=0  deg, 
<5=45  deg,  and  (p2  =  90  deg),  so  that  compression  in  Eq.  (24)  is 
normal  to  the  (110)  plane.  The  four  sets  of  initial  slip/twin  system 
resistances  listed  in  Table  7  are  investigated  for  each  orientation. 
For  comparison  purposes,  also  considered  in  what  follows  are 
predictions  for  a  fully  dense  polycrystal  (N=  100)  assuming  equal 
slip/twin  resistances. 

Axial  stresses  are  shown  in  Fig.  6(a).  Because  hydrostatic  pres¬ 
sures  are  nearly  identical  for  each  orientation  and  for  each  choice 
of  slip/twin  resistance  and  because  hydrostatic  pressures  are  rela¬ 
tively  large  compared  with  deviatoric  stress  components,  differ¬ 
ences  in  axial  stresses  among  various  cases  are  small,  at  most 
~5%.  Maximum  shear  stresses  computed  from  Eq.  (28)  are 
shown  in  Fig.  6(b).  Significant  differences  in  predicted  shear 
stresses  among  orientations  and  among  choices  of  initial  slip/twin 
strengths  are  evident.  Average  shear  stresses  supported  by  the 
polycrystal  fall  about  midway  between  those  for  the  stiffest  and 
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Fig.  6  Predicted  (a)  axial  stress  and  ( b )  maximum  shear  stress  for  AION  in  uniaxial  strain:  single 
crystals  (various  slip/twin  models)  and  a  fully  dense  polycrystal 


most  compliant  single  crystals  considered  in  Fig.  6(b).  For  single 
crystals  of  the  [110]  orientation,  lateral  stress  cr22^  o-33.  The  re¬ 
duction  in  slope  associated  with  pore  crushing  evident  in  Fig.  2 
(over  range  0.96>7>0.91)  is  absent  for  the  fully  dense  crystals 
considered  in  Fig.  6(a). 

Figure  7(a)  shows  the  predicted  twin  volume  fractions.  Clearly, 
evolution  of  twin  fractions  depends  strongly  on  initial  lattice  ori¬ 
entation  and  relative  hardness  of  slip/twin  systems.  For  the  single 
crystal  oriented  at  [100]  with  {111}  twinning  preferred,  the  twin 
volume  fraction  approaches  15%  at  7=0.9  (i.e.,  at  10%  compres¬ 
sive  strain).  For  a  single  crystal  with  the  same  set  of  slip/twin 
resistances  compressed  along  [110],  the  twin  volume  fraction  ap¬ 


proaches  only  3%  at  7=0.9.  Several  other  cases  including  single 
crystals  with  preferred  {110}  slip  demonstrate  almost  no  twinning 
(recall  also  that  twinning  is  prohibited  for  exclusive  {ill}  slip). 
Predicted  dislocation  densities  are  shown  in  Fig.  7(b).  Over  most 
of  the  compression  range,  the  [100]  orientation  with  preferred  slip 
on  {110}  or  {111}  planes  produce  the  largest  dislocation  density 
while  the  polycrystal  produces  the  lowest. 

Table  9  lists  for  7=0.97  and  7=0.9  (3%  and  10%  compressions, 
respectively),  the  ratio  of  maximum  shear  stress  on  each  of  the 
cleavage  planes  to  its  corresponding  theoretical  shear  strength 
given  in  Table  5.  Substantial  differences  are  evident  among  crystal 
orientations  and  initial  slip/twin  resistances  used  in  the  model.  For 


Fig.  7  Predicted  total  twin  volume  fraction  (curve  not  visible*-* negligible  values)  (a)  and  disloca¬ 
tion  density  ( b )  for  AION  in  uniaxial  strain:  single  crystals  (various  slip/twin  models)  and  a  fully 
dense  polycrystal 


Table  9  Ratio  of  maximum  shear  stress  to  theoretical  strengths  at  3%  and  10%  uniaxial 
compressions 


Crystal  orientation 

and  slip/twin  hardness  model 

Tmax 

/  Tf  {111} 

Tmax 

Itf  {110} 

rmax 

/tf  {100} 

7=0.97 

7=0.90 

7=0.97 

7=0.90 

7=0.97 

7=0.90 

[100]  all  systems  equal 

0.106 

0.119 

0.139 

0.157 

0 

0 

[100]  exclusive  {111}  slip 

0.099 

0.116 

0.130 

0.153 

0 

0 

[100]  preferred  {110}  slip 

0.063 

0.077 

0.084 

0.102 

0 

0 

[100]  preferred  {111}  twinning 

0.091 

0.105 

0.120 

0.139 

0 

0 

[110]  all  systems  equal 

0.116 

0.131 

0.140 

0.175 

0.234 

0.319 

[110]  exclusive  {111}  slip 

0.102 

0.125 

0.128 

0.224 

0.224 

0.435 

[110]  preferred  {110}  slip 

0.119 

0.142 

0.145 

0.191 

0.245 

0.346 

[110]  preferred  {111}  twinning 

0.093 

0.113 

0.121 

0.201 

0.217 

0.390 
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grains  of  the  [100]  orientation,  shear  stresses  on  {100}  planes  van¬ 
ish  identically.  Recovered  fragments  from  dynamic  experiments 
on  AlON  polycrystals  suggest  that  intragranular  fracture  occurs  on 
{111}  planes  [4],  At  7=0.9,  the  maximum  shear  stress  on  any  {111} 
plane  obtained  among  any  considered  set  of  orientations  and  ini¬ 
tial  slip/twin  resistances  is  ~14%  of  the  theoretical  strength.  Re¬ 
sults  in  Table  9  suggest  that  fracture  might  occur  first  on  {100} 
planes  for  crystals  of  [110]  orientation  if  the  theoretical  shear 
strength  computed  via  the  second  of  Eq.  (14)  is  realistic.  How¬ 
ever,  in  the  present  context,  Eq.  (14)  should  be  regarded  only  as 
an  order  of  magnitude  approximation  since  it  depends  only  on  the 
appropriate  shear  elastic  constant  and  interplanar  spacing  and  does 
not  account  for  details  of  atomic  bonding  (e.g.,  local  charge  im¬ 
balances  across  fracture  planes)  in  spinel  ceramics  such  as  AlON. 
The  experimentally  measured  fracture  toughness  for  MgAl204 
spinel  single  crystals  is  lowest  on  {100}  planes  [52],  often  quoted 
as  the  primary  observed  (tensile)  cleavage  planes  for  that  material; 
however,  atomic  modeling  [53]  has  suggested  that  {111}  planes  in 
spinel  ceramics  should  be  of  the  lowest  surface  energy  and  cleav¬ 
age  may  be  preferred  on  these  planes.  At  7=0.97,  the  approximate 
compressive  strain  in  an  AlON  polycrystal  at  its  HEL,  the  ratio  of 
maximum  shear  stress  to  theoretical  strength  on  any  of  the  planes 
considered  in  Table  9  is  substantial  (i.e.,  ~25%),  supporting  the 
notion  that  transgranular  fracture  of  grains  can  take  place  during 
high  rate  compression  [7,26], 

5  Conclusions 

A  model  has  been  developed  for  nonlinear  mechanics  of  elastic- 
plastic  ceramic  crystals  with  the  spinel  structure.  Nonlinear  elas¬ 
ticity,  anisotropy,  thermal  expansion,  dislocation  glide,  twinning, 
shear  fracture,  and  pore  collapse  are  addressed  for  AlON  crystals 
in  what  appears  to  be  the  first  documented  study  of  AlON  using 
large  deformation  crystal  plasticity  theory.  A  number  of  physically 
relevant  material  properties  corresponding  to  thermomechanical 
behavior  and  lattice  defects  in  AlON  have  been  compiled  or  com¬ 
puted.  Model  calculations  demonstrate  an  accurate  depiction  of 
pressures  and  shear  stresses  observed  experimentally  in  shock- 
loaded  polycrystals.  Key  findings  regarding  behaviors  of  AlON 
polycrystals  and  single  crystals  are  summarized  as  follows: 

•  Hydrostatic  pressures  and  axial  stresses  are  predicted  accu¬ 
rately  (the  latter  to  at  least  ~30  GPa)  when  nonlinear  elas¬ 
ticity  (i.e.,  pressure-  and  temperature-dependent  moduli)  and 
pore  collapse  are  incorporated  in  the  present  model.  Pres¬ 
sure  is  overestimated  when  pore  collapse  is  omitted  in  the 
model  and  axial  stress  is  underestimated  when  linear  elas¬ 
ticity  is  used. 

•  Because  of  the  large  number  of  potential  slip  and  twinning 
mechanisms  in  spinel  crystals:  slip  on  {111},  {110},  and 
{100}  planes  and  twinning  on  {111}  planes,  available  poly¬ 
crystal  test  data  are  insufficient  to  enable  a  unique  determi¬ 
nation  of  shear  strengths  for  all  inelastic  deformation  modes. 
Several  realistic  sets  of  initial  shear  resistances  have  been 
investigated,  as  has  a  model  for  shear-induced  fracture  dur¬ 
ing  impact  loading;  all  results  provide  a  reasonable  descrip¬ 
tion  of  macroscopic  axial-  and  shear  stress  data. 

•  Initial  resistances  to  continuum  slip  appear  comparable  to 
Peierls-Nabarro  stresses  computed  via  anisotropic  elasticity 
for  screw  dislocations  on  {111}  and  {110}  planes  in  the  rate- 
independent  limit.  Initial  resistance  to  twinning  on  {111} 
planes  appears  comparable  to  the  ratio  of  stacking  fault  en¬ 
ergy  to  the  magnitude  of  the  partial  Burgers  vector  associ¬ 
ated  with  twinning. 

•  Predictions  for  single  crystals  provide  estimates  of  effects  of 
lattice  orientation  on  the  response  to  loading  representative 
of  plate  impact.  Maximum  shear  stresses  attained  are  highly 
sensitive  to  crystallographic  orientation  and  to  initial  hard¬ 
ness  ratios  among  families  of  slip  and  twin  systems. 

•  Predicted  twin  volume  fractions  and  cleavage  fracture  ten- 

011013-10  /  Vol.  78,  JANUARY  2011 


dencies  vary  strongly  with  grain  orientation  and  initial  resis¬ 
tance  to  slip  or  twinning.  At  applied  strains  corresponding  to 
the  polycrystalline  HEL,  maximum  resolved  shear  stresses 
on  low-index  planes  can  approach  25%  of  the  theoretical 
strength. 

•  Results  point  to  the  need  for  experimental  measurements, 
either  in  situ  (since  AlON  is  transparent)  or  on  recovered 
samples,  of  twin  volume  fractions  to  quantify  preference,  or 
lack  thereof,  for  twinning  versus  slip  at  high  rates  of 
loading. 

The  content  of  the  present  paper  has  focused  on  constitutive 
model  development  and  predictive  calculations  for  AlON  single 
crystals  and  polycrystals.  Behavior  in  the  latter  has  been  described 
via  a  self-consistent  assumption  to  account  for  effects  of  porosity 
(Eq.  (21))  and  a  Taylor  assumption  for  deformation  of  crystals 
comprising  the  matrix  (Eq.  (24)).  A  more  refined  approach,  elimi¬ 
nating  the  above  assumptions,  would  fully  account  for  variations 
in  local  deformation  among  and  within  individual  grains  in  a  poly¬ 
crystal  as  well  as  possible  variations  in  behavior  of  individual 
pores  and  fracture  planes  in  the  microstructure.  Finite  element 
calculations,  whereby  individual  grains  and  grain  boundaries  are 
resolved  explicitly,  [15,44,45]  are  presently  in  progress;  results 
will  be  reported  in  future  work. 
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Appendix 

Expressions  for  shear  elastic  constant  A  used  in  expressions  for 
anisotropic  Peierls-Nabarro  stress  and  theoretical  strength  in  Eqs. 
(13)  and  (14)  and  Tables  4  and  5  are  listed  in  what  follows.  In  a 
coordinate  system  with  axes  aligned  along  the  cubic  axes,  the 
three  independent  second-order  elastic  constants  are  Cu,  C12,  and 
C44.  Components  of  the  second-order  compliance  tensor  are  [54] 

Hi  =  (C„  +  C12)/[(Cn  -  C12)(Cn  +  2C12)], 
i]2  =  -  Ci2/[(Cn  -  Ci2)(Cn  +  2C12)],  and  £44=1^44 

(Al) 

For  shearing  in  a  (110)  direction  on  a  {111}  plane,  the  appro¬ 
priate  shear  modulus  is  [55] 

A  =  l/{2[*44  +  (Hi  -  H2  -  2s44)/3]}  (A2) 

For  shearing  in  a  (110)  direction  on  a  {110}  plane,  the  appro¬ 
priate  shear  modulus  is  [54,56] 

A=l/[2(s11-Ji2)]  =  (C11-C12)/2  (A3) 

For  shearing  in  a  (110)  direction  on  a  {100}  plane,  the  appro¬ 
priate  shear  modulus  is  [54,56] 

/.r  =  i  /.s'44  =  C44  (A4) 
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